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The time evolution of spatial fluctuations in inhomogeneous d- dimensional biological systems is 
analyzed. A single species continuous growth model, in which the population disperses via diffusion 
and convection is considered. Time-independent environmental heterogeneities, such as a random 
distribution of nutrients or sunlight are modeled by quenched disorder in the growth rate. Lineariza- 
tion of this model of population dynamics shows that the fastest growing localized state dominates 
in a time proportional to a power of the logarithm of the system size. Using an analogy with a 
Schrodinger equation subject to a constant imaginary vector potential, we propose a delocalization 
transition for the steady state of the nonlinear problem at a critical convection threshold separating 
localized and extended states. In the limit of high convection velocity, the linearized growth problem 
in d dimensions exhibits singular scaling behavior described by a (d — l)-dimensional generalization 
of the noisy Burgers' equation, with universal singularities in the density of states associated with 
disorder averaged eigenvalues near the band edge in the complex plane. The Burgers mapping leads 
to unusual transverse spreading of convecting delocalized populations. 



PACS numbers; 05.70.Ln, 87.22. As, 05.40. +j 

I. LOCALIZATION AND POPULATION 
DYNAMICS 

The mathematical analysis of spatial patterns in bio- 
logical systems has been an object of intensive research 
for many years Both the dynamics and the equi- 

librium properties of certain model systems have been 
worked out in detail. Biological processes, such as the 
spread of a favored gene, population growth of species, 
ecological competition and so on, are often a combi- 
nation of diffusion and convection with some kind of 
back reaction. These systems are not conservative; both 
growth and death terms, possibly involving nonlineari- 
ties, change the number of individuals involved. A gen- 
eral form of such reaction-diffusion equations is M 



l^-f (v.V)c = f(c) 



DV^c, 



(1.1) 



where c(x, t) is the vector of reactants (e.g., species of 
bacteria, nutrients, etc.), D is a matrix of diffusivities, 
and f(c) describes the nonlinear reaction kinetics. The 
conservative term v • Vc represents a convective flux con- 
trolled by a drift velocity v, such as the flow of water in 
aqueous media, winds, etc. 

Although the literature discussing equations of this 
^> ; type is massive, the effect of spatial inhomogeneities in 
5^ the underlying medium is relatively unexplored A 
, disordered substrate may manifest itself in the above for- 
malism as quenched random diffusion constants, stochas- 
tic growth and death rates or randomness in the reaction 
term; it may reflect random concentration of environ- 
mental factors such as nutrients or toxins, or an inhomo- 
geneous illumination pattern projected onto, e.g., photo- 
synthetic bacteria. 

In this paper we study the effect of such heterogeneities 
in biological systems. As a model we take one of the sim- 
plest situations, the case of a single species, described by 



population number density c(x, t) , for which the reaction 
diffusion equation is a straightforward generalization of 
the Malthus-Verhulst growth model [0: 



9c(x, t) 
dt 



V • Vc(x, t) = DV^c{yi, t) + [a + C/(x)]c(x, t) 



-6c2(x, t), 



(1.2) 



where U{x) is a zero mean quenched random variable, 
and we take the convective velocity v to be constant in 
space and time. 

The homogeneous analog of this equation, without the 
convection term (i.e., Uhc) = and v = 0) was pro- 
posed by R. A. Fisher Q as a model for the spread of 
a favorable genetic mutation. It is also useful as a de- 
scription of a population dynamics, for which a is the 
difference between the linear birth and death rates, D 
reflects the effect of migration, and the term —bc^ repre- 
sents some self-limiting process, roughly proportional to 
the number of pairs of individuals at position x. In the 
non-random case ?7(x) = 0, there are two spatially homo- 
geneous fixed points: an unstable fixed point at c(x) = 0, 
in which there is no population at all, and stable fixed 
point at c*(x) = a/b, where the population saturates to 
the carrying capacity of the environment. Nonnegative 
initial configurations evolve smoothly toward the stable 
fixed point; analysis of the time development of spatial 
fluctuations in this model reveals that equilibrium can 
be reached via traveling solitonlike solutions, known as 
Fisher waves (l]§,§|. 

It is interesting to consider (|l.2| ) in the context of nu- 
cleation and spinodal decompos itio n. In the absence of 
convection, we can rewrite Eq. (|l.2|) as 



dc 
di 

where the potential is: 



-D 



.SF 
5^' 



(1.3) 
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FIG. 1. Schematic of a disordered biological substrate 
subject to diffusion and convection. The dark spots rep- 
resent environmental fluctuations, as exemplified by an 
homogeneous pattern of light projected onto a growing 
population of bacteria. 



F(c) = -(Vc)2 



$(c) 



with 



-J_[a + C/(x)lc2 + — 6c3. 
2D^ ^ 3D 



(1.4) 



(1.5) 



The time evolution of a small concentration inhomogene- 
ity thus seems to resemble the dynamics of an order pa- 
rameter quenched below its critical point and subject to 
a cubic Landau potential <i>(c) However, the usual 
additive Langevin thermal noise term is missing in the 
dynamics. This reflects the fact that our system admits 
an absorbing state, i.e., a population state which vanishes 
everywhere remains fixed at zero. A second crucial differ- 
ence arises because a smooth nonnegative initial condi- 
tion will remain nonnegative for all times |^ . The Fisher 
growth model is thus an unusual zero temperature "one- 
sided" spinodal decomposition problem. The one-sided 
nature of the dynamics insures that no unphysical effects 
arise due to the unbounded potential (1^) at large neg- 
ative c. Interesting studies exist of domain growth in, 
e.g., random Ising systems at zero temperature [||-@, 
but we are unaware of a similar body of work on "one- 
sided" population growth models. Another area of active 
research concerns zero temperature dynamical systems 
with multiplicative noise pd| . Here, however, the noise 
typically depends on both space and time, in contr ast to 
the purely space-dependent function C/(x) in Eq. (1.2). 



In Eq. (1.2), both convection and random fluctuations 
of the growth rate have been added to the original Fisher 
mod el. For the sake of concreteness, we can consider Eq. 
(1.2) as a model for a colony of bacteria which grows on 
an inhomogeneous substrate. The linear growth rate of 
the bacteria depends on a spatially random (but inex- 
haustible) food supply at each point, or other inhomoge- 
neous, time- independent, environmental factors such as 
the intensity of illumination. By reducing the light in- 
tensity, this linear growth rate could presumably assume 
both positive and negative values at different points in 
space. The food or illumination are quenched random 
variables. The bacteria diffuse, as well as undergo con- 
vection due to flow of the ambient medium with velocity 
V (Fig. 1) O. 



In this simple model of bacterial population growth, 
we neglect the dynamics of a diffusing nutrient supply or 
feedback from waste products ||l|. For growth in an in- 
homogeneous light source, for example, we might require 
a fixed homogeneous nutrient supply, possibly stabilized 
by large nutrient diffusion constants. Wakita et al. jl^] 
have studied the growth of Bacillus subtilis under various 
conditions, and found a large regime of low Agar density 
and high nutrient concentration w hich is well described 
by the homogeneous analogue of (1.2) without convec- 
tion, consistent with our assumptions. 



If the fluctuating growth rate ?7(x) is of order 5a, 
we shall see that convection will significantly perturb 
population growth when the flow rate v exceeds the 
corresponding change in the Fisher wave velocity 0| 
6vf — 2\/DSa. Upon taking as an effective diffusion 
constant for a motile bacteria D ^ 6 x 10~ cm /sec 
and growth rate fluctuations 5a ^ a = 10~^/sec, we find 
fiuctuating Fisher wave velocities of at most a few mi- 
crons per second. Thus, we shall be primarily concerned 
with small values of v and low Reynolds number fiows in 
our analysis. 



When randomness is introduced into (1.2), constant 
configurations are no longer possible steady-state solu- 
tions. Instead, the steady state is spatially modulated, 
reflecting the competition between the disorder and the 
diffusion term. Convection, on the other hand, tends to 
make the final state more uniform; For very high con- 
vection velocities, one might expect that each bacterium 
feels some "effective medium" average of the random en- 
vironment as it drifts rapidly from one site to the other. 
In this paper, we show that there is a phase transition 
from one regime to the other in the linearized growth 
problem and suggest that this sharp change of behavior 
persists in the steady state. 

In order to simplify our problem, let us assume that a 
stable nonzero steady-state population profile c*(x), ex- 
ists, as well as the unstable c(x) = steady state. The 
time evolution of small fluctuations around these (stable 
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or unsta ble) configurations is determined by the lineariz- 
ing Eq. (1^) near the fixed point functions. Linearizing 
about the "Gaussian" fixed point c(x) = leads to 

dc 

— ~L)V2c- V- Vc+ [a + C/(x)]c. (1.6) 

Linearization about the nontrivial stable steady-state 
configuration, c*(x), which satisfies 

DV^c* - V • Vc* + [a + f/(x)]c* - 2bc*^ = 0, (L7) 

leads to a similar equation for c'(x, i) = c(x, i) — c*(x), 
namely, 



'dt 



with 



ZJV^c'- V Vc' + [a + C/'(x)]c' (1.8) 



C/'(x) = C/(x) - 26c*(x). (1.9) 



We shall assume that c*(x) has no long range correla- 
tions, so that the linear operator in (1.9) now involves a 



redefined function C/'(x), with quenched random fiuctua- 
tions away from its mean value similar to those described 
byC/(x). 

One main interest in this paper is the evolution of 
small fluctuations about these fixed point configurations. 
Near the unstable c(x) = configuration, fluctuations 
grow until they reach the point for which the lineariza- 
tion fails. On the other hand, small perturbations of the 
stable state c*(x) will decay, so that the linear ap prox- 
i mat ion becomes better in time. Thus, Eqs. (L6) and 
(1.8), describe the short time growth of dilute popula- 
tions or the long time decay to the stable state. Although 
we shall frame most of our discussion in terms of the un- 



stable modes d escr ibed by (1.6), a very similar analysis 
applies to Eq. ( p^ ) . 

Similar equations emerge in a variety of physical situ- 
ations, such as models of chemical reactions and neural 
network [|l|. Miller and Wang et al. [|l5| have recently 
studied the spectrum of an operator describing diffu- 
sion of a passive scalar subject to a spatially random, 
time-independent velocity field, but without simple mul- 
tiplicative randomness. A closely related physical system 
from which we intend to draw, is fluctuating vortex lines 
in superconductors in the presence of columnar defects, 
with external magnetic field tilted away from the direc- 
tion parallel to the defects . The partition function 
Z(x, t) for a si ngle line at position x and height t then sat- 
isfies Eq. (L6) where D is given by temperature divided 
by the tilt modulus, D = T/{2ei), v is proportional to 
the tilt field, and U (x) corresponds to the columnar dis- 
order potential in the superconductor V^(x) normalized 
by the temperature, U{x.) = V{k)/T. Some basic facts 
about vortex lines are revived in Appendix A. It is inter- 
esting to note that Eq. (1.6) also describes the growth of 
monetary capital with diffusion, spatially varying interest 
rates and drift due to, say population migration. 



Another related problem concerns diffusion and drift 
of particles in a medium with randomly distributed traps 
]l8| pO[ . The long time decay of the density of active par- 
ticles as obtained experimentally from, e.g., photocon- 
duction studies in quasi one dimensional polymers p8| , 
is expected to exhibit stretched exponential relaxation in 
the absence of a bias 1 19 2^ . When a biasing electric field 
is present, the decay has a simple exponential prefactor 
pO| with however a transition from stretched exponential 
to exponential decay in the subleading behavior above a 
critical bias threshold |jl^. The coarse grained physics 
can be approximated by Eq. (1.2) without the nonlin- 
earity and U{x.) chosen so that all growth eigenvalues 
are negative when v = 0. Our primary concern here 
is with situations where at least some growth eigenval- 
ues are positive. However, a delocalization transition has 
also been invoked to describe an abrupt onset of a drift 
velocity as a function of the bias as the particle density 
decays in the trapping problem |^ . We discuss the rela- 
tion of the results of Ref. |l^ and ||2^ to the population 
biology problem treated here in Appendix C. 

After a Cole-Hopf transformation, i.e.. 



c(x, t) = exp[A$(x, i)] 



(1.10) 



the linear growth model described by Eq. (1.6) becomes 

(1.11) 



5t$(x, t) = L>V2$(x, t) + XD{y^f 
- V- V<I>(x, t)+a + U{x.), 



while c'(x, t) — ex p[A<l''(x, t)] generates a similar equa- 
tion from Eq. (|L8|) with J7(x) ^ U'{x). Chen et al. Q 
have proposed Eq. ( |l . 1 1| ) as a model for the dynamics of 
strongly driven charge-density waves wit h que nched dis- 
order. Later in this paper, we use Eq. (1.11) to study 
sample-to-sample fluctuations in /nc(x, t) in the limit of 
high convection velocities. We obtain exact analytic re- 
sults, which should be applicable both to linearized mod- 
els of population dynamics and to charge density waves. 
Equations (1.6) and (1 



1. 



may also be written as 



dtc = Cc, 
where the Liouville operator, e.g., 

/I^DV^-v- V + a + C/(x), 



(1.12) 



(1.13) 



generates the time evolution of the system. The spec- 
tra and eigenv alues of random non-Hermitian operators 
similar to Eq. (1.13) have attracted considerable interest 
recently |^^ . Provided linearization is an adequate ap- 
proximation, the dynamics of this system is determined 
by the eigenvalues and the eigenvectors of £. Near the 
stable flxed point c*(x) one expects only decaying modes, 
i.e., all real parts of the eigenvalue spectrum of £ are neg- 
ative, while near the unstable state c(x) = there are at 
least few positive, growing eigenstates. These expecta- 
tions can be demonstrated explicitly when randomness is 
absent, i.e., U{x) = 0. In this case the right eigenvectors 
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of the non-Hermitian operator ( 1.13| ) about the unstable 
fixed point are simple plane waves (?!>k(^) ~ e^^'^ which 
satisfy 

/:0«(x) = rk<^^(x), (1.14) 

with the complex eigenvalue spectrum 

T{k) ^ a-iv -k- Dk^. (1.15) 

The operator £ which corresponds to linearization about 
the nontrivial fixed point c*(x) ~ a/b has the same eigen- 
functions, with spectrum 



r'(k) 



iv-k- Dk' 



(1.16) 



Provided a > 0, the eigenvalues of (1.15) have a po sitive 
real part for small k, while all the eigenvalues of ( 1.16| ) 
have a negative real part. Note that the eigenfunctions 
are always delocalized plane waves. 

It is instructive to use the full nonlinear equation ( |1. 
(with C/(x) = 0) to trace the time evolution of a small 
random initial condition c(x, t = 0) in terms of the above 
eigenmodes when a > 0. Assume that initial conditions 
have contributions at all wave vectors k including k = 0. 
Although many eigenmodes are unstable near c(x) = 
and grow exponentially in time, the k = mode grows 
most rapidly. Once this uniform mode saturates at its 
fixed point value, it then acts to suppress all other grow- 
ing modes via mode coupling induced by the nonlinear 
mode coupling term —bc^(x.,t). These expectations are 
illustrated via a mean field solution of Eq. ( |l.2| ) in Ap- 
pendix B. 

The spatial characteristics of the eigenfunctions of the 
linearized growth operator change dramatically when 
randomness is present. When v = 0, the operator 



U (x) is Hermitian with real eigenvalues. 



and for strong enough disorder, all its eigenfunctions are 
real and localized; the localization length is smallest in 
the tails of the energy band, corresponding to extreme 
values of U{x.). In one and two dimensions, it is widely 
believed that all states are localized even for weak disor- 
der ||. 

When V ^ 0, the Liouville growth operator is no longer 
Hermitian, although it can still be diagonlized using a 
system of left and right eigenvectors. Since the convec- 
tion term in Eq. ( |l.(^ ) may be absorbed into the Lapla- 
cian by completing the square, V ^+21? ' right and 
left eigenfunctions of the new Liouville operator are re- 
lated to the eigenfunctions of £{v = 0) via an imaginary 
"gauge transformation" ; if 0„.v=o(x) is an eigenfunction 
of the Hermitian problem, then pM 



(/'n:v=o(a;) 



(1.17) 



are the eigenfunctions of the non-Hermitian operator 
with the same eigenvalue e„, up to a constant shift 20 1, 



en(v) = e„(v = 0) 



AD 



(1) 



provided that the localization length in the non- 
driven problem, is less then D/v. Thus, for small convec- 
tion velocities, there is a spectral rigidity - except for the 
shift, the real eigenvalue spectrum is locked to the values 
it had for v = 0. As v increases, however, some eigen- 
functions eventually become extended and the dynam- 
ics becomes sensitive to boundary effects. The correct 
eigenfunctions end eigenvalues are no longer related to 
the V = case by a simple transformation. In case of pe- 
riodic boundary conditions, complex eigenvalues and de- 
localized modes appear when D/v becomes smaller than 
As V is increased, these delocalized states appear first 
at the band center, for which the localization length is 
maximal, then move outwards. 

These expectations for the eigenvalue spectrum have 
been demonstrated by analytic work and a numerical 
analysis of a discrete lattice model, inspired by the 
physics of vortex lines ||2^. The corresponding lattice 
discretization of the nonlinear equation which motivates 
our present work reads 



^ = fEE[-"''''-x+e„(t) + e— c_„(t)] 

X V—1 

+ Y,[a + U{^)]c^{t)^hY,cl{t), (1.18) 



where Cx(i) is the species population at the sites {x} of a 
hypercubic lattice, and the are unit lattice vectors. 
Here w is proportional to the diffusion constant of the 
corresponding continuum model {w ^ D/£q), and g is 
proportional to the flow rate v. Up to subtraction pro- 
portional to I D in a, a, ?7(x) and h have the same in- 
terpretation as in the continuum model. When linearized 
about Cx = 0, Eq. (I.IS) may be written 



rfcx(t) 
dt 



5]£(x,x')cx'(i) 



(1.19) 



where the discrete Liouville operator C is the matrix 
d 



C 



-5][a + C/(x)]|x)(x| 



|x)(x + e,|] 



(1.20) 



As in the continuum case, the same equation arises when 
linearized about the a nontrivial fixed point c* provided 
we make the replacement U{k) — > U'{x) with ?7'(x) given 
by Eq. (|L9|) 

Typical spectra for a 1000-site model in one dimension 
with [/(x) uniformly distributed in the interval [—A, A] 
are shown in Fig. 2 for three values of g (x v [p5|. 
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FIG. 2. Energy spectra of one-dimensional 1000-site lat- 
tice model with randomness A/w = 1. The resulting 
spectrum for the same realization of the random poten- 
tial U{x) G [—A, A] is plotted here for three different 
values of g. (a) Case g < gi; all eigenstates are local- 
ized; (b) gi < g < g2', bubble of complex eigenvalues 
indicating extended states appears near the center of 
the band; (c) g2 < g; all the eigenstates are extended. 
(After Ref. [EHl). 



For g less than a critical value gi, all eigenmodes are 
localized, and the eigenvalues remain real and locked to 
their values for g = 0. For gi < g < g2, extended states 
with complex eigenvalues appear near the center of the 
band. Localized states still appear near the band edges. 
For 5 > (?2, every localized state is destroyed by the 
non-Hermitian perturbation, and all states are extended. 
In this limit, eigenfunctions are slightly perturbed Bloch 
states — the lattice version of plane waves. The spectrum 
is well approximated by the disorder-free limit, i.e., the 
lattice analogue of Eq. ( |1.15| ) 



T{k) = 2w cos(Ho -I- ig£o) 



(1.21) 



where £o is the lattice constant. 

With our definition of C, states near the top of the band 
should give a reasonable approximation to the spectrum 
of the continuous operator (1.14); Eq. ( 1.15| ) then de- 
scribes the upper edge of the ellipse of eigenvalues in Fig. 
2c. The states at the bottom of the band, on the other 
hand, have spatial characteristics which are artifacts of 
the lattice discretization. 

Figure 3 shows typical spectra for the discrete oper- 
ator in two dimensions Here, too, eigenfunctions 
near the top of the band should give good approxima- 
tion to growth modes in the continuum limit. Again, all 
eigenvalues remain real and the eigenfunctions remain lo- 
calized when g < gi. For gi < g < g2, however, extended 
and localized states coexist near the center of the band 
Ipsf . When g > g2, even most rapidly growing states at 
the top of the band are delocalized. For very large g (not 
shown) the spectrum resembles that of the disorder-free 
limit of the lattice model, similar to the one-dimensional 
case. However, as mentioned in Ref. p5| , this apparent 
simplification is actually a finite size artifact in d = 2: 
level repulsion leads to large modification of the Bloch 
wave functions and eigenvalues even for weak disorder in 
any sufficiently large system. Thus, chaotic spectra sim- 
ilar to Fig. 3c are always expected for large v due to 
disorder in the "thermodynamic limit" of large system 
sizes in 0? = 2. In Sec. 4 we argue that the growth for 
<? > <?2 when c? > 2 is in fact described by the {d — 1)- 
dimensional Schrodinger-like equation with space- and 
time-dependent randomness. The anomalous critical ex- 
ponents which describe this situation lead to universal 
power law singularities in the density of states near the 
band edge in the complex plane. 

If we take a gradient and set u (x, t) — — V(f>(x, i), 
the dynamical growth model ( |l.ll| ) which results from 
a Cole-Hopf transformation reads 

i9tu(x, t) + (v • V)u(x, t) + 2AL'(u(x, t) ■ V)u(x, t) 

= DV2u(x,i) -l-f(x), (1.22) 

with f (x) = VC/ (x) and subject to the constraint Vx u = 
0. This is a variant of the d-dimensional generalization 
of Burger's equation with noise studied 20 years ago by 
Forster et al. ||2^. In the form (1.11), such problems are 
sometimes referred to as "KPZ equations" . 
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FIG. 3. Typical spectra of the two-dimensional tight- 
binding non-Hermitian model with random site poten- 
tial, (a) Case g < gi\ delocalized states, real spectrum, 
(b) (71 < g < g2; extended states with complex eigen- 
values coexist with localized states with real eigenval- 
ues near the center of the band, (c) (j2 < fl; states 
at the tails of the band become extended, and there 
are complex eigenvalues even near the top of the band. 
(After Ref. [EHl.) 



In Sec. 5, we show that ( 1.22| ) is in the same universal- 
ity class as a simpler (d ~ l)-dimensional noisy Burgers' 
equation, and use this fact to study sample-to-sample 
fluctuations of Zn[c(x, t)] for large v. Thus, as the con- 
vective velocity v grows from small toward large values, 
the most rapidly growing modes evolve from a dynamics 
described by Anderson localization into a regime related 
to the Burgers' model for turbulence |p7[| . 

The remainder of this paper is organized as follows. In 
Sec. 2, we assume convective effects are small, and illus- 
trate the consequences of localization for "unbounded" 
population growth, i.e., growth at times before the non- 
linear term in Eq. (1.2) becomes important. In Sec. 3, 
we describe how these nonlinearities affect the growth 
when only a few modes are unstable relative to the state 
c(x) = 0. Unlike growth in a homogeneous environment, 
the most rapidly growing eigenfunction does not suppress 
all other unstable modes in this case. We suggest that the 
steady state undergoes a delocalization transition which 
occurs when the average growth rate or the convective 
velocity v are increased. 

In Sec. 4, we study the linearized growth problem 
in the large v limit, show that the physics is related to 
a (d — l)-dimensional Schrodinger equation and demon- 
strate that the average growth spectrum is singular for 
d > 2. This mapping leads to universal predictions for 
the randomness-dominated transverse wandering of a de- 
localized population as it drifts downstream. A d — 1- 
dimensional Burgers' equation is used to describe sample 
to sample fluctuations in /n[c(x, t)] in the limit of large 
V in Section 5. A number of related calculations are 
contained in three Appendices. Appendix C contains a 
discussion of the stretched exponential relaxation expect 
for populations in a medium with randomly distributed 
traps. 



II. UNBOUNDED GROWTH AND THE 
LOCALIZED LIMIT 

There is large literature on localization of electrons 
in disordered semiconductors |2^. Becaus e th e lin- 
earized growth modes which arise from Eq. (1.2) obey 



a Schrodinger-like equation, one might hope that results 
from electron localization theory would be applicable to 
the simple model of heterogeneous population dynamics 
discussed above. Of course, electronic states are either 
empty or full, according to Pauli exclusion principle. In 
contrast, many individual members of a species partici- 
pate in the growth modes determined by the continuum 
model equation studied here. Another differenc e is the 
form of the growth- limiting nonlinearity in Eq. (1.2). 



In this section, we explore the consequences of local- 
ization for the simple problem of "unbounded growth" , 
i.e., growth at times before the nonlinearity in Eq. ( |l.2| ) 
becomes important. We assume homogeneous initial con- 
ditions, weak convection, and parameters (small b, for ex- 
ample) such that the time domain over which unbounded 
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growth occurs is very large. Effects of the nonlinear term 
will be discussed in Sec. 3. 

Although many modes may now be growing exponen- 
tially, the fastest growing eigenfunction eventually dom- 
inates the center of mass of the evolving population dis- 
tribution. We show here that the time t* it takes for the 
"ground state" (i.e., the /asfesi growing eigenfunction) to 
win out in a large but finite domain grows very slowly as 
a function of the domain radius R. The precise form of 
t* (R) depends on the behavior of the density of localized 
states in the tail of the band of growth eigenvalues. For 
the simple discretized growth model discussed in Sec. 1, 
t*{R) ~ ln^^^(i?) where d is the dimensionality of space. 
Comparable results in semiconductors are usually deter- 
mined by electrons at Fermi energy in a partially filled 
band, and hence are less sensitive to the form of the den- 
sity of states. The time required for ground state domi- 
nance when populations grow in a spatially homogeneous 
environment is very different, t*{R) ^ R^. 

Giv en a time-independent Liouville operator C, as in 
Eq. ( 1.13 ), we describe growth in terms of a complete 
set of left and right eigenvectors {(/)^(x)} and {0^(x)}, 
with eigenvalues {r„} (the complex "energy spectrum"). 
The time evolution of c(x, t) is then given by 



</'n,v=o(x) 



-/^Ti |x — Xti I 



(2.4) 



c(x,i) = ^c„0,f (x) 



(2.1) 



where the coefficients {cn} are determined by the initial 
condition, 



Cn^ J d'^x (/)^(x)c(x,i 



0) 



(2.2) 



At long times the system will be dominated by the 
"ground state", i.e., the state for which the real part 
r„ is maximal. Throughout this paper, we assume that 
there is such a state, i.e., that the real part of the spec- 
trum is bounded above. The Liouville operator in Eq. 
(1.13) plays, with the replacement t —it and C —TL 
a role similar to the Hamiltonian TL in the Schrodinger 
equation ifidtc — TLc. We shall often use nomenclature 
from quantum mechanics, such as "energy spectrum" and 
"low-lying states" . However, because of the identification 
Ti. = — £, the ground state of the Hamiltonian is actually 
the state with the maximal eigenvalue of the Liouville op- 
erator, and the low-lying states of the Hamiltonian are 
those which grow fastest for the Liouville operator. 

We assume a small convective velocity, so that the left 
and right eigenfunctions are described by slightly dis- 
torted versions of the localized state for v = 0, as in Eq. 
(1.17), with normalization 



d-^x 0^,(x)0,^(x) 



(2.3) 



Deep in the band tail and close to the ground state, the 
localized eigenfunction for v will have the ap- 
proximate form 



where 6„ cx nt/"^ is a normalization constant and k„ is the 
inverse localization length associated with an eigenmode 
located at position x„. 

Of course, only the ground state is guaranteed to be 
strictly nonnegative pq |, as implied by Eq. (2.4). Or- 
thogonality with the ground state requires a small neg- 
ative part in the localized excited state eigenfunctions. 
Nevertheless, Eq. (2.4) should be a good approximation 




FIG. 4. The six fastest growing localized eigenfunctions 
for the two-dimensional lattice model (with v tx g = Q) 
which gave rise to the spectrum of Fig. 3(a). (Figure cour- 
tesy of Naomichi Hatano) 

for all modes in the tail of the growth spectrum, pro- 
vided the disorder is strong: Strong disorder means that 
all states are well localized, with approximately the same 
localization length. The six fastest growing modes on the 
square lattice which gave rise to the spectrum in Fig. 3a 
are shown in Fig. 4. 

Suppose the environment has been inoculated with 
a small uniform population of, say, bacteria at time 
t — Q. Then the projections {c„} of this initial condition 
onto the localized eigenmodes are approximately equal, 
c„ w Co for all n. The normalized species concentration, 



c(x,i) 



c(x,t) 



/ d'^x c(x, t) 



is given by 



c(x, t) 



^'f.(x)+E;>,?(x)e-(r--r")* 



(2.5) 



(2.6) 



where Ngs and the {Nn\ are normalization constants, 
TV,. = / d'^ 0« (x), 7V„ = / d'^ 0^(x). (2.7) 



We have separated out the "ground state," i.e., the most 
rapidly growing mode and the sums E„ are over the re- 
maining excited states. 
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The time necessary for the ground state to dominate 
the normaUzed species concentration in a large region of 



size R is clearly of order t* 



i/(r. 



i/Ar, 



where T* is the growth rate of the first excited state in 
the region. Let g{T) be the density of states per unit 
volume with growth rates in the interval between F and 
r + dV . Then, as we increase AF from zero, an excited 
state with gap AF will appear when 



i?''g(AF)AF w 1. 



(2.8) 



If (7(AF) were approximately constant near the band 
edge, then AF - l/R"^ and we would have t*{R) - i?''. 
However, there are in fact very few states in the tail of 
a band of localized eigenfunctions |^^. In Appendix C, 
we show that, for the simple lattice model of population 
growth discussed in Section 1, 



and 



+ A 



g(Ar) - exp[-(const./AF)'^/2] 



(2.9) 



(2.10) 



as AF 0. Solving Eq. ( |2.8| ) for large R now leads to 
AF ~ 1/ ln(i?/Zo)^^'' where Zo is the lattice constant of 
the model. The relaxation time for a system of size R is 
then 



t*{R)^W''{R/h) , 



localized growth . (2.11) 



The sparse population of growth rates near the ground 
state results in slow logarithmic growth of the relaxation 
time t*{R) with system size R. A very different size de- 
pendence results for the delocalized modes of the homo- 
geneous model. For delocalized plane wave eigenfunc- 
tions described by a spectrum like ( |1.15 ) , there are many 
more states close to the ground state. When v = the 
density of states 5(AF) - (AF)'*/^-! 

near the band edge 
in the delocalized limit, and the above argument leads to 



t*{R) ~ R^ 



delocalized growth. (2.12) 



III. INTERACTIONS AND A DELOCALIZATION 
TRANSITION 



IS 



A complete an alysis of the nonlinear "inte racti on" 
terms in Eq. (1.2), or its lattice equivalent Eq. (1.: 
beyond the scope of this paper. Some progress is possible, 
however, when only a few strongly localized growth eigen- 
values near the band edge of the Liouville operator have a 
positive real part. In this limit, it is easy to demonstrate 
that localized population dynamics with interactions dif- 
fers considerably from the dynamics of the plane waves 
which describe population growth in a homogeneous en- 
vironment. We shall also argue that the sharp mobility 
edge separating localized eigenfunctions from delocalized 



ones [ p5[ implies a delocalization transition with increas- 
ing convective velocity or average growth rate in popu- 
lation dynamics. This transition is clearly present in the 
linearized growth model, and we shall give arguments 
that it may be present as well in stea dy-st ate po pulation 
distributions described by Eqs. (F2) and ( l.lSj ). 
Let us write Eq. (1.2) in the form 



9c(x, t) 

at 



Cc{-K, t) - 6c2(x, t) 



(3.1) 



where C is the Liouville operator ( 1.13 ), and study the 
steady state which develops for long times. Upon ex- 
panding in the complete set of right eigenfunctions of £, 
with eigenvalues F„ 



c(x,t)=^c„(i)0^(x) 

n 

the dynamical equations read 
dcn{t) 



(3.2) 



dt 



where the mode couphng coefficients are 

Wn^mm'^b [ d-^X^^ (x)^« (x)0^, (x) 



(3.4) 



Upon combining the "gauge tran sformation" Eq. (1.17) 
with the approximate form ( ^.4| ) of th e lo calized eigen- 
functions, the orthogonality condition (p. 3D leads to p5[ 



I (2«„)'^ 

T{d)n{d) 



'J2«^ 

T{d)n{d) 



exp[v • (x - x„)/D - K„|x - x„| 



(3.5) 



exp[-v • (x - yin)/D - k„|x - x„ 



where 51^ is the surface area of a unit sphere in d dimen- 
sions. 

Now consider what happens in, say, the lattice pop- 
ulation growth model ( 1.18| ) when we vary the average 
growth rate a, starting with large negative values. As 
discussed in Appendix C, the real part of the growth 
spectrum broadens to a width of order wd + A about a 
due to diffusive hopping between sites and randomness. 
If —a ^ wd + A, then the real parts of all eigenvalues 
in Eq. (|3.3| ) are negative and the population becomes 
extinct. As we increase a however we eventually reach 
the simple but interesting situation where a small frac- 
tion of the eig envalues become positive. The coupled 
equations (B.2) then resemble the set of renormalization 
group recursion relations governing flows in the space of 
Hamiltonian coupling constants which describe equilib- 
rium critical points to one loop order p9| ]. Here, an ex- 
pansion in e = 4 — d is used to truncate an infinite set 
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of coupled differential equations. By analogy with criti- 
cal phenomena, we call the modes with positive Re(r„) 
"relevant variables" and those with negative Re(r„) "ir- 
relevant variables" . To a first approximation , we can 
simply discard the irrelevant variables in Eq. (3.3), be- 
cause their negative eigenvalues are not much affected by 
the dilute concentration of growing modes (see below), 
and they will eventually die out. 

The behavior of the relevant variables is also simple. 
Because these constitute only a small fraction of the total 
number of localized modes, they will be widely separated 
in space. The overlap integral which defines the coupling 
coefficients in Eq. ( |3.5| ) will then be negligible unless 
m = m! ~ n and the differential equations describing the 
localized eigenmodes decouple. 



dt 



r„c„(i) - WnC^it) 



where 



Wn = Wn,nn = b / d'^X^/.,^ (x) (x)]^ 



(3.6) 



(3.7) 



If the convective velocity is small, the eigenvalues r„ re- 
main locked at their values for v cx g = 0, and the ne- 
glected terms are smaller by a power of exp[— k/], where 
I is a typical spacing between relevant eigenmodes. 

The evolution of c(x, t) at long times (after the irrele- 
vant variables have died off) is determined by substitut- 
ing the solutions of Eq. (|3.6D , 



Cn{t) 



c„(0)e 



l + c„(0)f-(e 



r„t , 



1)' 



(3.8) 



into Eq. (3.2) and only summing over the unstable 
modes. In contrast to the homogeneous population dy- 
namics, where a single k = eigenfunction completely 
dominates the steady state (see Appendix B), the fastest 
growing localized eigenfunction does not interact appre- 
ciably with the other relevant variables. The steady-state 
fixed point c*(x) approached at long times is then char- 
acterized by many occupied modes. 



where the fixed point values c* are 



(3.9) 



(3.10) 



and Y^'^ means only unstable eigenfunctions make a 
nonzero contribution to the sum. This state is similar 
to the "Bose glass" phase of flux lines in Type II super- 
conductors, where a plays a role of a chemical potential 
for vortices Unlike the Bose glass, however, the 

number of degrees of freedom associated with each occu- 
pied localized state is highly variable. The "occup atio n 
number" Nn of an unstable eigenmode is, from Eq. (3.9), 



iV„ = < / d"^ 0^(x) 



(3.11) 



imr 




Imr 



Rer 



Rer 




i Imr 



Rer 



FIG. 5. Schematic of a one-dimensional eigenvalue spec- 
trum below and above the delocalization threshold. On 
the left, only a few localized states are unstable and the 
bubble of "irrelevant" delocalized states has little effect 
on the steady state. Spectra on the right have unstable 
delocalized states arising either from an increase in the av- 
erage growth rate a (top) or an increase in the convective 
velocity v with a — /AD held fixed (bottom). 

We checked the above analysis in one dimension for 
g (X V — by numerically determining the eigenfunctions 
in Eq. (1.20) for a particular realization of the random 
potential U{x) e [—A, A] with a sufficiently negative so 
that only four eigenvalues on a 700 site lattice were posi- 
tive. This situation is qualitatively similar to that shown 
on the left side of Fig. 5, except that there is no bubble 
of delocalized states in the center of the band for g = 0. 
However, provided Re(r„) is large and negative for these 
delocalized states, they rapidly die off and will not af- 
fect the fixed point describing the steady state. We then 
determined the steady stat e pop ulation distribution for 
the full nonlinear equation (1.18) under the same condi- 
tions. This steady state is compared with the four rele- 
vant eigenmodes of the linearized problem in Fig. 6. The 
peaks in the exact steady state (top part of the figure) 
do indeed occur precisely at the locations of the four un- 
stable growth eigenfunctions. It is easy to see from Eq. 

(U) that Wn,mi — const. Kn ■ Sincc K„ (the inverse 
width of the localized eigenfunctions) is approximately 
independent of n and equal to the lattice constant in the 
t ail o f the band (see, e.g., Ref. Q), it follows from Eq. 
( |3.10| ) that the heights of the peaks in the steady state 
are proportional to the growth eigenvalues of the unstable 
modes in this simple model of population dynamics. We 
have checked that this relationship between peak heights 
and growth rates is satisfied by the steady state of 
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steady State 



0.0 



0.0 



350.0 



700.0 



0.50 



First Four Eigenstates 








j.t \ 


1 i 



0.0 350.0 700.0 

FIG. 6. (Top) Steady-state population distribution for 
a 700-site tight-binding model with site random growth 
rates and nonlinear interactions. (Bottom) Plots of the 
four unstable eigenfunctions obtained by linearizing the 
Liouville operator for the same realization of the random 
potential about the state of zero population. The posi- 
tions of these eigenfunctions match perfectly the peaks in 
the steady state of the nonlinear problem. The heights 
of the peaks in the steady state are proportional to the 
eigenvalues of the corresponding eigenfunctions. 

Fig. 6 with an accuracy of a few perc ent. Thu s, a pop- 
ulation of bacteria described by (|l.2|) or (I.IS) evolves 



toward a steady-state distribution given by the ground 
state and first few excited states of a Schrodinger-like 
equation. 

We next discuss the case shown schematically on the 
right side of Fig. 5, where a number of delocalized states 
have become unstable. This regime could be accessed 
either by increasing the mean growth rate a, or by in- 
creasing V which enlarges the bubble of extended states. 
(To pin the center of the band one should actually in- 
crease V while holding a — /AD fixed). Because the 
Liouville operator is real, these delocalized modes occur 
in pairs: if (f>n{^) is a mode with eigenvalue r„, then 
(f>^{x)* is an eigenfunction with eigenvalue T* . As shown 



in detail by Brauwer et al. , these modes are (at least 
for weak randomness) approximately plane waves, char- 
acterized by nonzero wave vector pairs k and —k. The 
dynamics changes dramatically as soon as the first pair of 
delocalized eigenfunctions becomes relevant. Now, there 
will be nontrivial mode couplings between the newly un- 
stable delocalized modes and each other, as well as with 
the unstable localized ones discussed earlier. By analogy 
with the physics of tilted vortex lines interacting with 
columnar defects , we now expect macroscopic oc- 

cupation of modes near q ~ 0, similar to Bose-Einstein 
condensation. A broad background extended species pop- 
ulation should now be superimposed on the peaks which 
represent localized modes, as indicated schematically in 
Fig. 7. The proportion of, say, bacteria incorporated into 
this background should increase with increasing a. De- 
localization arises because populations can drift between 
growth "hot spots" instead of dying out when the con- 
vection velocity is high. According to the mode couplings 
in the Fourier basis displayed in Eq. (B.3) of Appendix 
B, fixed point value for the q = mode is determined by 
the values of the unstable pairs (k, —k). 

There should be a large difference in the response of 
the steady state to a small change in v for localized and 
delocalized steady states. The spectrum "unzips" fur- 
ther into the complex plane with increasing v oc g as 
more modes delocalize. For spectra like that on the left 
side of Fig. 5, these additional delocalized modes are ir- 
relevant and should not affect the steady state. The only 
change occurs due to the dist orti on of the relevant local- 
ized modes according to Eq. (|3.5| ) . For spectra like those 
on the right side of Fig. 5, however, increasing v (with 
a — /AD constant) leads to more relevant delocalized 
modes, with large changes in the corresponding steady 
state. Note that the "carrying capacity" defined by Eq. 
11) diverges when g k„ and a mode described by 



(I 



an eigenfunction like (3^) become s d elocalized. How- 
ever, coupling coefficients such as (3.7) remain finite at 
the delocalization transition. 

Although we have discussed the dynamics using one- 
dimensional spectra, we expect similar delocalization 
phenomena with increasing v or a with two-dimensional 
spectra such as that in Fig. 3b. If a is adjusted so that 
only a few localized modes are relevant, and the growth 
modes are strongly locaHzed, the steady state should look 
like Fig. 4 with peak heights proportional to the growth 
eigenvalues. With increasing growth rate or convection 
velocity, eventually both delocalized and new localized 
modes should start to participate in the steady state. 
The coexistence of localized and extended eigenmodes at 
the same real part of the energy in d = 2 was discussed in 
p5[ : anisotropic localized eigenfunctions for g = will 
delocalize sooner if their most extended direction coin- 
cides with the direction of g. We expect that the result- 
ing extended state wave functions are streaked out in the 
direction of g, and that these streaks will be reflected 
into the steady-state species population. 
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a 




FIG. 7. Schematic "phase diagram" indicating regimes of 
locahzed steady states (lower left) and extended steady 
states (upper right) as a function of the average growth 
rate a and the convection velocity v. (A region of large 
negative a where all initial populations eventually become 
extinct is not shown.) 

It is possible that a delocalization transition arises even 
for V = with increasing a (or increasing diffusion con- 
stant D). As the number or spatial extent of the relevant 
eigenmodes goes up, eventually the localized states be- 
gin to overlap. The steady state may become "extended" 
in this limit, in the sense that individual members of a 
species can hop easily from one growth hot spot to an- 
other. Again, the response of the steady state to a small 
change in convection velocity (away from zero) could be 
an indicator of this transition. A sharp transition is sug- 
gested by analogy with the "Bose glass" transition which 
occurs for increasing magnetic fields for vortices in Type 
II superconductors with columnar defects (see Appendix 
A and Ref. 0). 

A more detailed discussion of these interesting delocal- 
ization transitions will be presented in a future publica- 
tion |32|. However, in the next two sections, we shall 
make some progress in describing the low-lying states 
which describe the growth modes in rf = 2 for large con- 
vective velocities. 



we mean velocities so large that even the "ground state," 
i.e., the growth mode with the largest eigenvalue, is de- 
localized. 



(4.1) 



where is the localization length of the most rapidly 
growing eigenfunction. In this regime, it is helpful to 
exploit the analogy with the equilibrium statistical me- 
chanics of a 2 -|- 1-dimensional vortex line described in 
Appendix A. In this analogy, a vortex trajectory rep- 
resents the path taken by a particular growing popula- 
tion of, say, bacteria. As is evident from Fig. 10, when 
ti ^ oo the component of the "magnetic field" tilting 
the equivalent elastic vortex line leads to configurations 
more nearly perpendicular to the columnar defects repre- 
senting the disordered growth rates. The transverse fluc- 
tuations of the tilted vortex line's trajectory are unim- 
peded in the r directions, while in the xy plane it sees the 
cross sections of the columnar defects. One might guess 
that in this large tilt limit the vortex simply wanders 
diffusively along the r direction, but acts like a directed 
polymer in a 1 -I- 1-dimensional medium with point-like 
disorder when projected into the xy plane. More gener- 
ally, the physics of tilted vortex lines with columnar pins 
in d -|- 1 dimensions should be related to directed poly- 
mers in (d — 1) -|- 1 dimensions with point-like disorder. 
A great deal is known about these problems |3^ . 

To apply similar ideas to linearized population growth 
models, we assum e that v is in the x direction, and 
rewrite Eq. (1.16) as 



dtc + vdxC - 



Ddlc 



Dy\c+[a + U{x,ri_)]c (4.2) 



where rj^ represents all spatial coordinates perpendicu- 
lar to X. Specifically, we explore the long-time dynamics 
generated by this equation with a delta function initial 
condition. 



lirn c{x, r_L, t) 



5{x)5''-^{v^) 



(4.3) 



corresponding to a point inoculation of population at the 
origin. We look for a delocalized solution valid for long 
times in the limit of large v. The overall exponential 
time dependence generated by a, the constant part of the 
growth rate, and the diffusion with drift we expect in the 
(x, t) variables may be incorporated via the substitution 



IV. BURGERS' EQUATION AND THE LIMIT OF 
LARGE V 



A. Qualitative Discussion 

We now study "unbounded growth" in the limit of 
large v, with the goal of better understanding com- 
plicated spectra like those for the two-dimensional lin- 
earized growth problem shown in Fig. 3c. By large v 



c{x,Y^_,t) 



e W{x,Ti) (4.4) 



where W{x, y^_) is to be determined. Note that c(x, rj^, t) 
becomes proportional to 5{x) as t — s- 0, so the initial con- 
dition ([4.3[) requires 



lim Vl^(.T,r_L) = <5''"^(r_L) 



(4.5) 



Upon inserting Eq. (4.4) in ([l.2|), we find 
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DdlW{x, r^) + DV\W{x, r^) 
+U{x,Y^)W{x,ri_) . (4.6) 



According to the ansatz (4.4), c{x,v^,t) is only appre- 
ciable for 



vt\ < 2VDt 



(4.7) 



so the second term on the left-hand side of (4^) is smaller 
than the first by a factor of order jtv^^ and can be ne- 
glected in the limit of long times. In the remaining equa- 
tion for W(x^ r^), which has no explicit time dependence, 
we expect that the term d'^W{x, r±) can be neglected for 
large x and t compared to the single x derivative which 
appears on the left-hand side. The resulting equation is 
an imaginary time Schrodinger equation, where x plays 
the role of "time" , 

vd,Wix,r^) « DVlW{x,r^) + U{x,r^)Wix,r^) . 

(4.8) 

Note that the random "potential" U{x,Yx) depends both 
on the "time" x and on the additional (d — 1) spatial di- 
rections. 

Consider the application of this mapping to two- 
dimensional species populations with strong convection. 
Assume for simplicity that a > 0, so that the popula- 
tion grows on average as it convects and diffuses down- 
stream. For fixed x, the solution W{x,y) of the result- 
ing (1 -f l)-dimensional S chro dinger equation (subject to 
the boundary condition (45)) describes the distribution 
in 2/ of a growing species population which has traveled 
through random distribution of growth rates for a time 
of order t = x/v. The results of extensive studies of 
the Schrodinger equation in 1 + 1 dimensions with a 
space and time-dependent random potential may be in- 
terpreted as follows: For any fixed y value, imagine trac- 
ing the genealogy of, say, all bacteria which have reached 
particular position {x,y). As x — *■ oo, the overwhelm- 
ing majority of bacteria near the point (x, y) will have 
evolved along a spatially convoluted optimal path of espe- 
cially favorable growth rates. The fraction of bacteria 
whose ancestors come along this route is favored over all 
routes in W{x,y) by a exponential factor exp[c'a;'^], 
where c' is a constant. The exponent lo (which describes 
the fluctuations in the ground state energy in the anal- 
ogous problem in the statistical mechanics of flux lines) 
in known to be w = 1/3 exactly |Q. Any particular 
path y{x) of optimal evolution itself wanders with typi- 
cal transverse fluctuations which behave like 



X 



C 



(4.9) 



where ( = = 2/3 is a universal critical exponent, 
independent of the exact probability distribution of the 
randomness and other details. Because 2/3 > 1/2, this 



optimal path will be well defined even in the presence of 
diffusion. A schematic of a set of optimal paths is shown 
in Fig. 8a. The population distribution c{x,y,t — x/v) 
is sketched in Fig. 8b. The exponent C, also controls 
the overall transverse spread in y of the spatially varying 
population for fixed x. This "superdiffusive" spreading 
(C > 1/2) arises because trajectories which would be rare 
in conventional diffusion lead to strong amplification if 
they pass through regions of particularly favored growth. 
The exponent uj determines the size of the fluctuations 
in c{x,y,t = x/v), which ride on top of an overall expo- 
nential growth, c{x, y,t — x/v) ~ e' 



j: / V 



Similar behavior is expected for the 2 -I- 1-dimensional 
random Schrodinger equation which results for convect- 
ing populations with randomness in d = 3, with the uni- 
versal exponents C w 0.59 w 3/5 and lu = 2( — 1 = 1/5 
p3| . The growing population again becomes streaked-out 
in streamwise direction, but with a nontrivial wandering 
transverse to the stream. 




In c(x, y, t = x/v) 




(b) 

FIG. 8. (a) Trajectories for a growing species which will 
produce a particularly large population at the point {x, y) 
at time of order t = x/v. The population associated with 
a given point is spread over a region with typical size 
Dx/v. For large x, bacteria which have traveled along 
such a path of favorable growth rates will dominate the 
population at {x,y). (b) Schematic of the logarithm of the 
population c{x,y,t — x/v) discussed above as a function 
of y. Typical fluctuations in population size away from 
simple exponential growth are of order exp[c'a;"], and the 
population from an initial source of point inoculation has 
spread out a distance of order x'' . 
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B. Average Growth Spectrum 

We now study large v growth spectra like those in Figs. 
2c and 3c averaged over many realizations of the disorder. 
Reference presents numerical evidence and qualita- 
tive arguments that the effect o f diso rder on the spectrum 
of the non-Hermitian operator ( 1.13|) in one dimension for 
large v is in fact very weak |3j ] . It was argued that the 
more complicated chaotic eigenvalue spectra observed in 
d = 2 were due to level repulsion of discrete eigenval- 
ues in the complex plane. Here we first show explicitly 
that perturbation theory in the disorder simp ly le ads to 
a "free particle" dispersion relation like Eq. (LIE) with 
renormalized values of a, v, and D in one dimension. 
We then demonstrate that the same perturbation analy- 
sis is singular in higher dimensions, consistent with the 
mapping onto a (d— l)-dimensional random Schrodinger 
equation described above. 

For a given configuration of growth rates, we apply 
standa rd sec ond-order perturbation theory |^5[ to the op- 
erator (1.13). The resulting growth rate spectrum, start- 
ing with a plane wave set of basis functions reads. 



e(k)=eo(k) + (k|C/|k)+^ 



l(q|t/|k)p 

eo(k) - eo(q) 



0[C/3(x)] 
(4.10) 



With V in the x direction, we avera ge over disorder and 
take the limit — > oo in Eq. ( 4.10| ), and find 



e-(fc,,fci) = eo(fcx,fc±) + l^i^' J ^ 



Jd-l 



1 



i2n)d-Uv{q, - k,) + D{q^ - k^) 

(4.17) 

where Iq is the lattice constant of a lattice and w e hav e 
kept only the small k expansion displayed in Eq. ( 4.12| ). 

Evaluation of the average spectrum is particularly sim- 
ple in one dimension, because the transverse wave vectors 
kj_ and are absent. Upon extending the integration 
limits on to ±oo, setting qx — kx = Px and symmetriz- 
ing in Px , we have 

e{kx) = a - ivkx - Dkl + ^k^^ 



dpx 



D 



2tt D^pI ~ {iv + 2DkxY 
loAH 



Dki 



6{iv + 2Dkx) 



(4.18) 



Expanding the disorder correction in k leads to renor- 
malized values of a, w, and D 



where for a basis appropriate to the lattice growth model 
(1.20) with TV sites, we have 



-y 



e*''-''|x), 



and 



(4.11) 



VR 



Dr = D + 21qA'^D^/3v^ 



(4.19a) 



(4.19b) 



(4.19c) 



with 



eoik) = a~iv -k- Dk^ + 0{k^) , (4.12) 



We see that disorder increases the mean growth rate, de- 
creases the effective convective velocity, and increases the 
effective diffusion constant in an expansion in the dimen- 
sionless ratio IqA^D/v^. 

The new growth eigenfunctions are slightly perturbed 
plane waves [p5[ 



The averages over disorder are easily calculated by first 
using this lattice model, and then passing to the con- 
tinuum limit. With C/(x) uniformly distributed in the 
interval [—A, A], we have 



1 



C/(x)f/(x') = -A2^,^, 



(k|C/|k) = 



and 



|(q|C/|k)|2 = A737V 



(4.14) 



where the overbar represents a disorder average. It fol- 
lows that 



(4.15) 



(4.16) 



\ky^\k) + J2 



{q\U\k) 



q^k 



(4.20) 



and are clearly still delocalized. Equation ( 4.191 ) sug- 
gests that localization sets in whenever 



> 1, 



(4.21) 



i.e., localization occurs when the effective drift velocity is 
renormalized to zero. The physical basis for this criterion 
is as follows. In the absence of convection, the maximum 
growth rate associated with the operator (1.13) (i.e., the 
ground state eigenfunction of the corresponding Hamilto- 
nian) will be given by minimizing the variational function 
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W[c{x)] 



I 



H 



dc{x) 
dx 



\u{xy{x) 



dx 



J c'^{x) dx 



(4.22) 



The random potential U (x) in the continuum hmit now 
has correlator 



1 



Uix)U{x') = ^A'^loS{x ~ x') 



(4.23) 



where Iq is a microscopic cutoff, of order the l attice con- 
stant in a discrete growth model like (1.20). We as- 
sume an exponentially localized nodeless "ground state" 
growth eigenfunction centered on the orig in cq (x) 
exp(— K|a:|), and replace the random part of ( 4.22|) by its 
root-mean square value. Upon neglecting dimensionlcss 
coefficients of order unity, we find 



Win) - Dk^ - A(/ok)^/^ 
which is minimized for k = Kq? with 

.o«(A/yv^)'/^ 



(4.24) 



(4.25) 



and we again neglect constants of order unity. The 
"gauge transformation" (1.17) allows this state to remain 



localized only if vjD < kq, which is equivalent to the 
criterion ( 4.21 ). 

We now demonstrate the singularities which arise for 
large v||x in higher dimensions d>2. Upon setting 



Px = qx 



(4.26) 



and symmetrizing in px , Eq. ( 4.17 ) becomes 



where a is a positive dimensionlcss coefficient. Thus a 
"free particle" spectrum of the form 

eo{kx, ky) « an - iVRk^ - D^kl - D^kl (4.30) 

is only a good approximation provided v is large enough 
such that 



-^L„ < 1 



(4.31) 



Equation (4.31) is consistent with numerical results and 
a criterion based on a level repulsion argument for the 
lattice model However, for any fixed value of v, 

there will always be nontrivial changes in the growth 
spectrum for sufficiently large Ly, consistent with chaotic 
spectra like that exhibited in Fig. 3c. The mapping onto 
the physics of a (d — l)-dimensional random Schrodinger 
equation suggests that, when evaluated to all orders 
in perturbation theory, the renormalized wave vector- 
dependent transverse diffusion constant actually diverges 
as (7j_ ^ HI), D^{qj_) ^ q]^'' ^, so that the disorder 
averaged renormalized spectrum takes the form, valid for 
small wave vectors. 



With C = 2/3, we have 
e{kx,ky) an- ivakx - 



A\k 



i/C 



(4.32) 



^x^x ^ A±ky^^ 

(two dimensions) . (4.33) 



e(fcx,kj 



Dxki 



D<kj 



a — ivkx - 
dpx f d'^^^p^ 
2tt J (27r)''-i 
D,pl + D^{pl-kl 



[D^pI + D^{p\ - kl)Y - pl{iv + 2Dxkxf 



(4.27) 



where we have separated the diffusion term in eo(fcx, fc_L) 
into components parallel and perpendicular to v with 
Dx = D± = D. It is tedious but straightforward to 
demonstrate that a, v, and suffer only finite renor- 
malizations similar to Eqs. (4.19) when this formula is 
expanded in kx- However, upon setting kx = and we 
find a diverging renormalization of D±^ in the long wave- 
length part of the integral. 



0^1 = 0^ 



12v 



(27r)'i-i pI 



less singular terms. 

(4.28) 



Consider the meaning of this infrared divergence for 
d = 2. In a finite system of spatial extent Ly <ti Lx in 
the y direction, we have 



1 



A2/2 

a L, 

vDjL ^ 



(4.29) 



This growth spectrum implies that length scales in the 
transverse direction scale with the 2/3 power of stream- 
wise length scales, consistent with Eq. (|4.9| ). The expo- 
nent 3/2 controlling transverse fluctuations also appears 
in the renormalization group treatment of fluctuations 
in lnc(r, t) presented in Sec. 5. In d = 3, the mapping 
onto a random Schrodinger equation leads to an exponent 
C ~ 3/5, and the small wave vector form, 

e{kx,k±) w a_R - ivRkx - D^k^ - A±k^^ 

(three dimensions). (4.34) 

It is interesting to com pute the consequences of spec- 
tra like ( 4.33| ) and ( 4.34 ) for a density of states in the 
complex plane, deflned by 



5(ei,e2) 



dkx 1 ^ .-J^ j-r 

^ J J2^^^'' 
Re e{kx,k±^)]S[e2 - Im l{kx,k^)]. 



(4.35) 



Note that this is the density of states associated with 
disordered averaged eigenvalues. If the fluctuations of 
the eigenvalues away from their average values are small, 
this quantity will be the same as the disordered average 
density of states studied, e.g., in Ref. p2|. We measure 
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energies relative to an so that the top of the band (cor- 
responding to the ground state of the equivalent Hamil- 
tonian) occurs near the origin. Straightward calculations 
then lead to the predictions 

gi^i,<^2) = 0, ei > or el > -eivji/D^; 

otherwise (two dimensions) (4.36) 



and 



9(^1,^2) = 0, 61 > or ej > -eivj^/D^; 

g(6i,62)«(|6i|-Df6>2)l/5^ 

otherwise (three dimensions) (4.37) 

close to the band edge. 

These results are valid for small real and imaginary 
energies 61 and 62. The density of states thus vanishes 
outside a parabolic boundary, as in Fig. 3c and similar to 
its behavior in a pure system. However, it diverges as the 
boundary is approached in d = 2 and vanishes in d = 3 
with new universal critical exponents when disorder is 
present. Disorder thus has a strong influence in the ther- 
modynamic limit for large v when d > 2, even though all 
states are delocalized. In a homogeneous growth model, 
this density of states has a square root divergence near 
the boundary in d = 2 and approaches a constant in 
d = 3. 

As a crude approximation, one could use the average 
growth spectrum to estimate the behavior of, for example 
c(a;, y, i) in d = 2, 

c{x,y,t) sa 



'y_pe{kx,ky)t pik^x+ikyy 

2tt I 2tt 



e^-^'-y^e— • -""cik^, ky, t ^ 0) (4.38) 



where c{kxTky^t = 0) is the Fourier transformed initial 
condition. However, it is easier and more systematic to 
study the statistics of ln[c(a;, t)], as is often case for sys- 
tems with multiplicative noise. This is done in the next 
Section. 



V. SPACE-TIME FLUCTUATIONS OF Zn[c(x,t)] 

In this section we give a detailed analysis of the re- 
sponse of a homogeneous biological system to the in- 
troduction of quenched random inhomogeneities in the 
growth rate. Let us look, for concretcness, at a homo- 
geneous initial condition, i.e., c(x, = 0) = const., which 
then evolves under the influence of some kind of quenched 
spatial randomness. As discussed in Section I, if the ini- 
tial, constant density c(x, i = 0) is very small compared 
to a/6, the short time dynamics of this system is deter- 
mined by Eq. (|1.6D, and for c(x, t) « c*(x), i.e., near the 
stable fixed point, the long-time decay into c* is given 



by Eq. {^^- The dynamic renormalization group ap- 
proach presented here is thus directly relevant to the de- 
cay of small unstable fluctuations into the steady state 
or to the growth of unstable modes if b is small enough 
such that the long-time behavior of the linearized prob- 
lem inanifests itself before the nonlinear reaction term 
in ( |l.2| ) becomes important. Moreover, as discussed in 
Section III, the basic features of the linearized problem, 
such as its eigenvalues and eigenstates, can be used in 
certain limits to assess the nonlinear time evolution, as 



inEq. (^. 

We assume the term —bc^ in is negligible, and 

use the fact that c(x, t) is always nonnegative to define 
its logarithm $(x, t), via the transformation 



c(x,i) =e5iJ*(>'^*)+'^*c(x,t = 0) 



(5.1) 



With v||x, the function $(x, t) then satisfies (for a uni- 
form initial species population) 



5t$ -I- vdx<^> = DV'^<S> + ■^(V$)2 + C/(x) , 



(5.2) 



with the initial condition ^(x, t = 0) = 0. An identi- 
cal equation was studied numerically and via a scaling 
ansatz by Chen et al. as a model of charge density 
waves interacting with quenched disorder. Here, we use 
the renormalization group to study this problem analyt- 
ically using the methods of Ref . . 

We have singled out the direction parallel to the drift 
as a;, while the d — 1-dimensional perpendicular space 
will be denoted r_L, as in Sec. 4. Equation ( ^.2|) takes the 
form 

dt^{Y^,x,t) + wc):r«'(r±,x,i) = 

i^iVi$(rj_, X, t) + D,dl^{v^,x, t) (5.3) 
+ Ax(Vi$)2 + K{d^<^f + U{v^,x) , 

with = Dx = D, and X± — = A. The random 
function U{r±,x) satisfies 



U{r^,x)Uir'^,x) = T5{x ~ a;')*^"*"^ [r_L - r^] (5.4) 

where the correlator strength is related to the spread 
[—A, A] of growth rates in a lattice model by T oc A^£q. 
We now impose a change of scale 



r_L srj. 

X s^x 
t sH 



(5.5a) 
(5.5b) 
(5.5c) 
(5.5d) 



where s is a renormalization group scale factor. 

Under this scale transformation the parameters of Eq. 
(^.3D change according to 
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s2 



(5.6) 



V 



V 



If the nonlinearities are absent, i.e., At 



A, 



0, Eq. 



( p. 3]) becomes exactly solvable. In this case, "I>(rj_,a::,t) 
is given by 



, du) f°° dK 



(5.7) 



gikj_ ■rj_ gifc^Sg-i!^* 



where the Fourier transform of $(rj^, i), satisfies 

$(ki,fcx,t^) = G'o(k_L,fcx,w)t/(k_L,fcx,t^) (5.8) 

Note that since the disorder is time independent, its 
Fourier transform satisfies 

?^(q±,9x,w) = ?7(q,9x)27r(5(w), 

and 

((7(k, fc,)[/(k', fc^)) = T(27r)'^<5(ki - k'^)5(fc, - fci). 

(5.9) 

The bare propagator Go(kj^, fc^, tj) is 

1 



Go(kj_, k^,,uj) 



(5.10) 



In this case the nai ve s cahng treatment gives exact re- 
sults. /,From Eqs. (5.6) one finds that keeping D±, A 
and V fixed under the scale transformation requires 



2, 



r, = 2, 



3-d 



(5.11) 



dimensions to calculate the average. It is straightfor- 
ward to show that the simple scaling of the linear theory 
sketched above implies that W{Lx, L±) takes the form, 
for d = 2 



W{L,,L^) = L^hiL^/Li) 



(5.12b) 



where x = 1/2 and C = 1/2. To obtain (|5.12bD , u se Eg . 
( ^.23 ) below with the Gaussian exponents frora (fe.lll) 
The same results follow directly from Eq. (5^) Plf 



qa:<L 



dqx 
-1 27r 



dq± 



(5.12c) 



which also leads to the conclusion that h{x) ^ x for large 
X. Chen et al. simulated the full nonlinear equation (5.2) 
in d = 2 and found x = 0.5 ± 0.05 and C = 0.85 ± 0.05. 
Below we determine these exponents exactly. 

To consider systematically the effect of the nonlinear- 
ities, we perform a perturbative expansion around the 
exact solution embodied in (5.7- 5j-0| ) in powers of A_l 
and Xx- We begin by rewriting ( ^.7^ as an integral over 
the momentum up to some cutoff: 



$(rj^,x,t) 



dij 
2^ 



27r 



(a) 



(5.13) 



The term proportional to D^. then scales to zero, con- 
firming that it is negligible at long wavelengths and low 
frequencies compared to the single derivative in the con- 
vective term. The same analysis suggests that the non- 
linear coupling Xx is also irrelevant, A^; s~^^+'^^/^Aa;, 
while A_L — > s^~''/^A_L so that d = 3 is the upper critical 
dimension for this problem. 

These results become easier to understand if we con- 
sider their impact on a concrete physical quantity. Fol- 
lowing Rcf. pl| , we consider the sample-to-sample fiuc- 
tuations in ^Jr±,x,t) in different random environments 
all with the same physical dimensions Lx and L±, 



W{Lx,L^)^<P^x,r^,t) 



ln^[c(x,r±,t)e-"*] . 



(5.12a) 



For population dynamics of, say, bacteria, one could di- 
vide a single large colony into many patches with these 



(b) 



-o 



y y = — 



<-|k(k-q)-^-k(k,-q,) 



FIG. 9. ( a) Dia grammatic representation of the integral 
equation (5.14). (b) Iteration solution of Eq. (5.14) and 
this meaning of the vertex in this diagrammatic series. 
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FIG. 10. (a) One-loop correction to the self-energy ob- 
tained by averaging over the noise, (b) Perturbative cor- 
rection to the two-point function <l?(x, t)^{0, 0). (c) Three 
diagrams which contribute to the one-loop effective vertex 
renormalization. 

Here, A is a cutoff initially of order i^^, where is a mi- 
crosc opic length scale of order the lattice constant. Equa- 
tion (5.3) now becomes 



U{q±,qx,uj) 



-oo 27r 



dqx 



27r 7o (27r)'^-i 
-q_L ■ [k± - q]_L 
$(k_L - q_L,fc^ -qx,uj- 9)<^{ci, qx,^)] 



-^q_L ■ [k± - q]_L - ^qx[kx ~ fe] 



which can be represented gr aphi cally as in Fig. 9a. The 
formal iterative solution of ( |5.14| ) is shown in Fig. 9b. 

The self-energy S(kj_, /c^;, w) (Fig. 10a) yields the 
renormalized velocity vr and diffusion constants and 
D^. The renormalization of T is represented by the 
graph in Fig. 10b, and vertex renormalization graphs 
for \x and Aj_ are shown in Fig. 10c. The corrections 
to A_L, Aa:, Dx and v do not diverge in the infrared limit 
in any dimension; naive perturbation theory yields, how- 
ever, infrared divergent corrections to D and T for d < 3 
as one takes the limit k — > 0, 



dqq'^'^ (5.15) 



5T = Kd-i 



f A^T 

\8D^\v\ 



dqq 



d-4 



(5.16) 



where Kd ~ Sd/i^n)'^ and Sd is the sur face area of the 
d-dimensional unit sphere. Equation ( 5.15|) is similar 
to the result (4.28) obtai ned b y cons iderin g the average 
growth rate. Equations ( ^.15 ) and (5.16) confirm that 
the upper critical dimension of theory is d = 3, below 
which renormalization group techniques are needed to 
take care for the infrared divergences of the loop inte- 
grals. Upon carrying out the procedure of Ref. ||2^], we 
find that the renormalization group flow equations for 
the relevant variables take the form 



dD±_ 
dl 
dA 
Hi 

dX± 

dl 

dv 

m 



3-d 



Di_ z-2 + Kd-i , , , 
A 2z~d+l-2a-'q 
A^ [a + z - 2] 
v[z — rj] 



Kd- 



(5.17) 
(5.18) 
(5.19) 
(5.20) 



We have set s = e~^, reduced the cutoff from A to e~^A, 
and the corrections to the naive scaling results ( ^.6D are 
proportional to the dimensionless coupling constant. 
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AA^ 

2m^\ 



(5.21) 



The couplings Dx and A^^ again scale to zero, even at the 
nontrivial fixed point discussed below. 

We now set rj = z and a = 2 — z to ensure that A_l and 
V remain unchanged by our renormalization procedure. 
Using Eqs. ( 5.17)-(5.20) we can calculate the flow of the 
coupling constant g, 



dg 
'dl 



Kd 



2d -5 
4d-4 



9 ■ 



(5.22) 



(5.14) The flxed point g* is obtained by taking ^ = 0. For 
d — 2 (i.e., two-dimensional disordered growth model 
with one parallel and one perpendicular direction) the 
Gaussian fixed point at = is unstable, while g* = 
^ ^ ^ is an attractive fixed point which corresponds to 



a — 1/2, z ~ 3/2, and t] = 3/2 (ry should be equal to z 
since there are no infrared diverging correction to v) . 

We now integrate the recursion relations until £ is large, 
so that all transients have died off. The renormalization 
group homogeneity relation for the mean square fiuctu- 
ations in ln[x{z,r±,t)e~°'*] in a system with dimension 
Lx and Lj_ then reads 



WiLx,L^)=e 



W{Lxe 



(5.23) 



and 



The prefactor follows from Eq. (5.5d), with r/ 
a = 2 — z; we expect these results to be correct to all 
orders in perturbation theory due to the Galilean invari- 
ance of the problem after averaging over disorder [p3| . 
Upon choosing I — I* such that LxC'^^* = 1, Eq. (5.23|) 



can be rewritten in a scaling form similar to (5.121) 



where 



W{Lx,Ly)^Ltf-^hiL^/Li) 



C=l/z 



(5.24) 



(5.25) 



We have thus confirmed the scaling ansatz Eq. ( 5.5b| ) 
PH with the specific predictions r/ = 1/z = 2/3 and 
X — 4^ — 2 — 2/3. These exact exponents differ some- 
what, however, from the numerical estimates of Chen et 
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al. it would be interesting to see if the agreement 
improves with larger system sizes. 

For d = 3, there is no perturbatively acce ssible fixed 
point of the one loop recursion relation ( |3.22|) . However, 
extensive numerical work exists suggesting a stable non- 
trivial fixed point with z w 5/3. Thu s we predict from 
Eq. ( 5.23| ) that the scaling relation ( 5.5b| ) holds, with 
( ^ 1/z = 3/5 and X = 4C - 2 = 2/5 in d = 3. 

Our analysis of the evolution of ln[c{x, t)] shows in ef- 
fect that the asymptotic behavior of the nonlinear equa- 
tion (5.2) has the same critical exponents as a conven- 
tional noisy Burgers equation [23] with one less dimen- 
sion. The same dimensional reduction for critical expo- 
nents d — > {d— 1) was described for the original linearized 
growth model in Sec. IVa. A related dimensional reduc- 
tion has been found by Tang et al. |3^ for a model of 
driven depinning in anisotropic media, and by Obuhkov 
for directed percolation JSTj . 

It is interesting to comment about the role of v in this 
perturbative renormalization calculation. The small di- 
mensionless parameter of the series, g, is proportional to 
1/v and diverges as u — s- 0. When w — > the theory is in 
the strong coupling limit where we expect localized states 
in the band tail. A signal of this phase transition is the 
(finite) one loop correction to v obtained from Eq. (5.14) 



Sv 



dqq 



d-2 



(5.26) 



Although this correction is finite it suggests the existence 
of a critical value of v below which vr is zero, consistent 
with localized states unaffected by convection. 
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APPENDIX A: POPULATION DYNAMICS AND 
VORTEX CONFIGURATIONS 

In this appendix we review the statistical mechanics 
of a superconducting vortex line with a columnar pin- 
ning potential p7| , and show that its partition function 
evolves with sample thickness in the same way as the lin- 
earized population dynamics problem studied here. The 



different configurations of the vortex line, described here 
simply as an elastic string, are related to possible space- 
time trajectories of populations which diffuse, grow and 
drift in an inhomogeneous but time- independent environ- 
ment. 

Consider a superconductor sample of thickness L, 
pierced by "columnar pins," which are long aligned 
columns of damaged material, illustrated schematically 
in Fig. 11. The vortex can usually be described by a sin- 
gle valued trajectory r(r), where we assume defects are 
aligned with the r direction, r = x x y. The free energy 
of this problem may be written as: 




FIG. 11. Vortex line in a superconductor with columnar 
disorder. If H_i_ = 0, the flux line is localized, i.e., trapped 
by one or few pins into some region in the plane perpen- 
dicular to the correlated disorder. The coordinate r along 
the columns plays role of time. As the external magnetic 
field is tilted away from the columns, the flux line tends 
to delocalize and tilt in the direction of the external field. 



F[rir)] 



2 JO 
1 
2 



dr 



dr(r) 
dr 



dr V[r{T)] ~ n 



VP 
'At: 



dr 



(A.l) 

dr ' 

dz 



where e is the tilt modulus of the flux line and the elas- 

2 



tic contribution 



is the first nontrivial term in 



the small tipping angle expansion of the line energy of a 
nearly straight vortex line. V{r) is the random potential 
which arises from a r-independent set of disorder-induced 
columnar pinning potentials (with its average value sub- 
tracted off), and Hj_ is a perpendicular magnetic field. 

The partition function 2^(x, r; xq, 0) associated with a 
vortex which start at position xq at the bottom of the 
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sample (r = 0) and terminates at position x at tempera- 
ture T somewhere in the interior at height r is given by 
a path integral: 



Z(yi, r;xo,0) 



■(r)=x 



I?r(r)l 



-^[r(T)]/r 



(A.2) 



r(0)=xo 



Standard path integral techniques pq ] may be used to 
show that Z(x, r; xq, 0) obeys the Schrodinger-like equa- 
tion, 

- T— = V^Z - -hi • WZ + -^Z + V{x)Z 
OT 2e e 2e 

(A.3) 

where hj_ = is the dimensionless perpendicular 

field. Thus, the growth of the partition function Z(x, t) 
of a flexible line with height r is t he same as linearized 
population growth model Eq. ( |l.6| ), with the identifica- 
tions T -^t, D T/e, V -> h±/e, U{x) V{x)/T and 
a^hl/2eT. 
Because 



lim Z(x, r; xo, 0) = (5 (x — xq) 

T— >o 



(A.4) 



the full partition function Z( x, t;xo,0) is in fact the 
Green's function for Eq. (1.6), assuming a delta func- 
tion initial condition of population at position Xq and 
time i = r = 0. 



dcajt) 
dt 



acQ 



q#0 



(B.5) 



q#0.q'#k 



Cq(t)ck-q' (t) 



(B.6) 



The mean field approximation consists of neglecting 
terms in th e mo des with k 7^ 0, shown in brackets in Eqs. 
(B.5) and ( |B.6| ). The approximate differential equations 
which remain have solutions 



co(i) 



Ck(i) 



[1 + Mii(eat_i)] 



Ck(0)e 



[l + i|i(eat_i)p 



(B.7) 



(B., 



where cq = Oco(t — 0). The time evolution of c(x, i) 
is now co mple tely determined by inserting these results 
into Eq. (B.2). Note that there are many growing modes 
with k 7^ at short times when a > 0. However, when 
t ^ ln(6co/a) the denominators of Eqs. (B.7) and ( p3.8| ) 
cause co(i) to approach cq = a/b, and all modes with 
k 7^ to decay away with the spectrum (1.16). 



APPENDIX B: MEAN FIELD THEORY OF 
HOMOGENEOUS POPULATION DYNAMICS 



We start with the homogeneous analogue of Eq. (1.2) 



d 

— c(x, i) + V Vc(x, t) 



(B.l) 



D V^c(x, t) + ac{x, t) - bc^ (x, t) 
Upon decomposing c(x, t) into Fourier modes 



(B.2) 



where $7 is the volume of a c?-dimensional box with peri- 
odic boundary conditions, we have 

"^'"^^^ = TkCk - ^ ^Cq(t)Cq,(i)^k,q+q' (B.3) 



dt 



q,q' 



with the complex growth spectrum 



a + «v • k - Dk^ 



(B.4) 



In the spirit of Bogoliubov approximation for the collec- 
tive excitations in superfluids |3^ , we separate out the 
k = mode and write 



APPENDIX C: DENSITY OF STATES NEAR 
THE BAND EDGE 

In this appendix we calculate the density of states 
(DOS) in the tail of the band of growth rates, i.e., near 
the ground state. We set the growth bias a = 0, since it 
does not affect the statistics of the DOS. The discussion 
below assumes localized eigenmodes in the "relevant" or 
unstable part of the spectrum, where the convection term 
g cx V just produces a trivial shift in the eigenvalues, and 
hence we neglect g as well. 

Consider, then, a d-dimensional hypercubic lattice 
with edge length Iq, where the potential energy at each 
site is in the range C/(x) g [—A, A]. Unbounded proba- 
bility measures may give different results. 

The continuum Liouville operator approximated by the 
lattice model may be written as: 



C = DV"^ + L/(x). 



(C.l) 



It is easy to see that the DOS function p{E) in this model 
is bounded from above by A, such that p(E) — > as 
E ^ A. The tail of the DOS is determined by the 
range of energies in which the DOS is determined by 
rare events, characterized by large spatial regions with 
low potential energy [M . 



19 



Let us estimate these fluctuations in the following way: 
the probability to find a hypersphere of radius R which 
contains only blocks of potential energy U larger than 
A — Vo is approximately 



P{R, U> A~Vo) 



2 A 



(R/ta) 



exp 



In 



2A 



(C.2) 



The energy of a state confined by this rare fluctuation is 
given approximately by 



EkA~Vo- D/R^, 



(C.3) 



so that the probability to get an energy between E and 



E + dE using a sphere of radius R is p{R, E) ~ -gy^ | Vo=-Ei 



I.e. 



p{R, E) ~ exp 



Ry , (\ E 

— ln\ 

h) V2 2A 



2i?2A / 



(C.4) 



This expression is well defined for 



< i? < OO. 



Optimizing p(i?, E^ with respect to R gives, up to log- 
arithmic corrections, a maximum at the lower limit 
R* « \JD/{A - E), so that as £; A from below, 
p{E) = p{R*, E) vanishes according to 



P^E)^eM-iD/il[A~E]) 



d/2i 



(C.5) 



Now the DOS is proportional to p{E), i.e., at the tail 
of the distribution we have g{E) ~ gQp(E), where 50 is 
some normalization of order the DOS in the middle of 
the band. 

Let us consider now the tight-binding analog of the 
above model. The on site potential is U{x), taken from a 
square distrib ution in the range [—A, A]. The Liouvillian 
is, from Eq. (|l.20|). 



+ ^C/(x)|x)(x| . 



(C.6) 



The eigenenergies of this Hamiltonian are bounded, — A— 
w < En < A + w with w ^ D/£q. The states in the 
tails correspond to rare spatial fluctuations of U (x) . The 
probability to find such fluctuations (e.g., a region of ra- 
dius R in which the on site potential is within a specified 
energy interval Vb of the maximum value A) is the same 
as in the previous model. The energy spectrum of such 
fluctuation is given approximately by 



A - Vo + w 2J cos{k^£o), 



(C.7) 



where k 
relation 



1/R. Thus, states in this tail obey the 



A + w-Va-w/R^ 



(C.8) 



so that a result of the form (C.5) is applicable here also, 
with the energy measured from the edge of the band de- 
fined by the lattice model ||3^ . 

We can use the density of states result (C.5) and the 
eigenvalue spectra displayed in Figs. 2 to understand 
the results of Refs. |]l9| and j2^ for particles diffusing 
and drifting with random traps in one dimension. We 
start with the expansion of c(x, t) in the complete set of 
eigenf unctions of Eq. (1.2) with 6 = 0, as in Eq. (2.1). 
For V = g = 0, all eigenvalues r„ are real and negative 
for this problem. As discussed in |20|, the eigenvalues 
close to zero arise from the rare regions discussed above 
for the density of states. The spectrum for the lattice 
model will look like Fig. (2a), with however the top of 
the band just touching the origin. Upon assuming a uni- 
form initial condition c(x, 0), we take all Cn ~ const and 
integrate over space to obtain 



d'^xc{x,t) (X / g(r)e"dr 



(C.9) 



Where g{r) is the density of states. At i — > 00, the be- 
havior of the total number of surviving particles is dom- 
inated by energies near the top of the band, where 



5(r)cxexp[-|ro/r| 



<i/2l 



(C.IO) 



A saddle point evaluation of Eq. (C.9) in d dimensions 
then leads to 



7Vtot(i) ocexp[-(t/to)^], 



(C.ll) 



in agreement with Refs. ||19| and ||20|] . 

For V = g < gi, the effect of nonzero convection is sim- 
ply a rigid downward shift of the spectrum, by an amount 
of /AD 1^^. When g > gi in one dimension, the bub- 
ble of delocalized states shown in fig. (2b) will appear in 
the center of the band of negative eigenvalues. However, 
the behavior of the density of states at the top of band 
is unchanged, and we find 



V t 

Ntot{t) oc exp[-— ]exp[-(t/to) 



1/31 



(C.12) 



in agreement with Ref. ||T^. Delocalization affects the 
long time decay only when g > g2, when all states are 
delocalized, as in Fig. 2c. The density of states in one 
dimension is then 



g{AT) cc 1/{AT) 



1/2 



(C.13) 



which leads an additional exponential contribution to the 
decay of Eq. (C.9) (up to logarithmic corrections), in 
qualitative agreement with the transition as a function 
of the bias found by Movaghar et. al. m^. It would be 
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interesting to use the spectra displayed in Figs. 2 and 
3 to extract the short and intermediate time behavior of 



as well as the effect of drift in higher dimensions. 



e.g. V. Privman, Tr ends in Stat. Phys. 1 (1994); J. Cardy,^ 
cond-mat /96071 63 ; J. Cardy and Uwe C. Tauber, cond 



mat/9609151 



Note that the delocalization transition which describes 
particles diffusing and converting in the presence of traps 
occurs at the top of the band. For population biology 
problems, one must consider positive growth eigenval- 
ues and the phenomena of interest typically occur for 
gi < g < g2- The localization transition of interest to us 
in this paper occurs when the mobility edge crosses the 
origin, as in Fig. 5. 
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